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£C) ■ We propose a dispersive representation of the charged pion vector form 

f**"*- ' factor that is consistent with chiral symmetry and fulfills the constraints 

. imposed by analyticity and unitarity. Unknown parameters are fitted to 

the very precise data on r — > ty ir°v T decays obtained by Belle, leading 
to a good description of the corresponding spectral function up to a nn 
squared invariant mass s ~ 1.5 GeV 2 . We discuss the effect of isospin 
CO ' corrections, and obtain the values of low energy observables. For larger 

values of s, this representation is complemented with a phenomenological 
description to allow its implementation in the new TAUOLA hadronic 
currents. 
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1 Introduction 



Last years have witnessed a notorious progress in the knowledge of the two-pion 
system, both from theoretical and experimental sides. In particular, the high 
precision measurements of the charged and neutral pion vector form factors per- 
formed at the flavour factories BaBar [TJ, Belle 0, CMD-2 [3], KLOE @] and 
SND [5] have significantly improved the accuracy of previous data. As it is well 
known, theoretical predictions for these form factors cannot be obtained analyt- 
ically from first principles through standard calculations, since this involves in 
general the hadronization of QCD currents in a nonperturbative energy regime 
(E < 1 GeV). In order to overcome this problem, one can rely on effective mod- 
els [6j [7J [8] that intend to describe the involved dynamics within a simplified 
scheme. In general, in these models the form factors are required to satisfy 
the proper behaviour at very low energies, and the effect of intermediate reso- 
nances is taken into account through the inclusion of Breit-Wigner functions. 
In this way one can obtain phenomenologically adequate hadronic matrix ele- 
ments [TUl [TTJ [12] , which have been included e.g. in the standard Monte Carlo 
Generator for tau decays, TAUOLA [13] ■ However, in general these models in- 
clude several ad-hoc assumptions, and can even have problems of consistency 
with QCD [141 115 ) . In addition, even when the agreement with experimental 
data can be very good, usually the model parameters can be hardly related to 
the underlying strong interaction theory. 

Another possible approach is to consider just the symmetry properties of 
QCD in order to build a general effective action adequate for the desired energy 
regime. At very low energies E <C M p , where M p is the /o(770) resonance mass, 
the approximate chiral symmetry of QCD allows to build the effective quantum 
held theory known as Chiral Perturbation Theory (x?T) [15]. The latter pro- 
vides a successful description of the low-energy phenomenology of strong and 
electroweak interactions, in which hadronic observables are calculated through 
an expansion in powers of ratios of momenta and masses of the lightest degrees 
of freedom (light pseudoscalar mesons) over a chiral symmetry breaking scale, 
AttF^ ~ 1.2 GeV. However, for E ~ M p the expansion parameters become large, 
and new degrees of freedom, namely the lowest-lying light-flavoured resonances, 
become active. Even though in this regime there is no straightforward expan- 
sion parameter, one can build an effective theory by considering an expansion 
in powers of the inverse of the number of colours, 1/Nc, with the introduction 
of resonances as active fields in the effective action. Indeed, it is found that 
this approach allows to describe satisfactorily most salient features of meson 
phenomenology [17l [18], which suggests that the large- Nq limit of QCD is a 
good starting point to derive a chiral Lagrangian that includes resonance fields 

pi na eu. 

In this work we study one of the simplest hadronic observables, namely the 
pion vector form factor Fy(s), defined through 

<7r°7r-|dV U |0) = V2F^(s) (p n - - p T of , (1) 

where s = q 2 = (p n - +p 7T o) 2 . For s > 0, this form factor is probed by the decay 
t~ — > Tr~Tr°is T , while in the isospin symmetry limit it can be experimentally 
measured from e + e~ — » tt + -k~ (s > 0) and elastic e~7r + scattering (s < 0). 
The analysis of Fy(s) allows to increase our knowledge of the hadronization of 
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QCD currents in the intermediate energy region, where the presence of meson 
resonances plays a crucial role. On the other hand, the analysis of isospin 
breaking corrections to the form factors in tau decays and e + e~ scattering [22, 
|2"31 [Ml |2"5] is essential for the evaluation of the hadronic contribution to the 
anomalous magnetic moment of the muon a^, which provides a stringent test 
of new physics |26j . 

The theoretical analysis of the pion vector form factor has been addressed by 
several authors in the last years. At very low energies, Fy(s) has been computed 
in xPT up to 0(p 6 ) [37J[2H]- Then, to enlarge the domain of applicability up to 
~ 1 GeV, unitarization techniques [29] [30] and dispersion relations have been 
employed [3T] [32] [33]. Moreover, in order to go beyond this energy region, 
the inclusion of the p(1450) resonance [31] and even a tower of resonances, 
inspired in the Nc — > oo limit [35] [3B], have been proposed. Our work is a 
sort of extension of those in Refs. [3T] [33] [37] [35], in which the authors analyze 
■kit and Kit vector form factors considering 0(p 4 ) expressions obtained from a 
chiral effective theory that includes the dominant resonance exchange, followed 
by an Omnes-like resummation of final state interactions. Our procedure is 
similar to that proposed in Ref. [35] for the Kw vector form factor: we consider 
an n-subtracted dispersion relation for Fy(s) in which the input elastic phase 
shift 5{(s) is taken from the effective theory, resumming the chiral loops into 
the denominator of the 0(p 4 ) form factor. This ensures to fulfill unitarity 
and analyticity constraints. It is seen that a phenomenologically good result 
is obtained with three subtractions, hence our expression for Fy(s) depends 
on four parameters: M p , F^, and two subtraction constants ot\ and a>2 (the 
remaining subtraction constant is fixed by the normalization of the form factor). 
These constants can be related to chiral low-energy observables [27] [28] , namely 
the squared charged pion radius (r 2 ) v and the coefficients of 0(s 2 ) and 0(s 3 ) 
terms in the chiral expansion, c v and d v , respectively. 

The above described approach is able to provide a good description of the 
r~ — > n~n®v T spectral function for a squared irir invariant mass up to about 
Smax — 1.5 GeV 2 . Beyond these energies, we propose a complementary expres- 
sion for the form factor that includes the effects of the excited states p' and 
p" , matching smoothly the previous one at s ~ s max . Our result for the full 
form factor can be useful to improve the new version of TAUOLA 39 , which 
presently includes the expressions obtained within a chiral Lagrangian frame- 
work with resonances gD] HU S3] [H] . This is important [15] not only for 
the proper simulation of backgrounds and subsequent signal extraction at the 
more frequent tau decay modes but also for the analysis of rare processes and 
the searches of new physics [15] . 

The article is organized as follows: in Section [5] we obtain a dispersive rep- 
resentation of Fy{s) and discuss the inclusion of isospin breaking corrections. 
The model parameters are fitted to experimental data up to s ~ 1.5 GeV 2 , 
and our input for 5\{s) in the elastic region is confronted with present experi- 
mental values. Then, in Section [3J we extend our parametrization of Fy{s) to 
higher energies, including the effective contribution of excited resonances. The 
agreement with experimental data is shown. In Section [U we present the results 
for the low-energy observables related to our subtraction constants. Finally, in 
Section [5] we state the conclusions of our analysis. 
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2 Low energy description of Fy(s) 

As stated, at very low energies the pion vector form factor is well described by 
xPT. Let us first consider the limit of exact isospin symmetry. At 0(p A ), one 
has 07J 



F v (s) 



F2 



96n 2 F 2 



(2) 



where Lg(p) is one of the renormalized low-energy coupling constants in the 
chiral Lagrangian. The functions Ap(s,p 2 ) are given by 



A P (s, p 2 ) = log ^ + 8^ - £ + o*(e) log 
p z s 3 

where the phase space function ap(s) reads 



CTp(s) + 1 
tJp(s) - 1 



(3) 



Ms) = V 1 - 4 — 



(4) 



On the other hand, the computation of Fy(s) from a chiral Lagrangian that in- 
cludes the lowest-lying vector meson multiplet as active resonance fields yields [19] 



F v (s) 



1 



Fy Gy 
F 2 



My 



(5) 



where My — M p , and Fy and Gy measure the strength of the pV p and pirir 
couplings, respectively, being the quark vector current. This tree level result 
corresponds to the leading term in powers of 1/Nc, and it is 0{p A ) in the chiral 
expansion. If the form factor is required to vanish in the limit of large s, then 
one gets the relation FyGy = F 2 , which yields 



F V (S) 



M 2 P 
Mj-s 



(6) 



Comparing with Eq. ([2]), the low energy xPT coupling Lg is predicted to be 



La = 



FyGy 

2M 2 



F 2 
2Mj 



7.2- 10" 3 



(7) 



in very good agreement with the value obtained from phenomenology. This 
shows explicitly that the p(770) contribution is the dominant physical effect in 
the vector form factor of the pion. Now, as stated in Ref. [3T] , one can do better 
and match Eq. ^ to the 0(p i ) xPT result in Eq. @, including the final state 
interactions encoded in the chiral loop functions Ap(s,p 2 ): 



F V (S) 



M 2 



M 2 



96n 2 F 2 



Ms) + \Ak{s) 



(8) 



We omit from now on the explicit dependence on the p scale, taking p = M p . 
The results do not depend significantly on changes in this scale. From Eq. ©, 
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unitarity and analyticity constraints lead to the Omnes exponentiation of the 
full loop function [3"T] . 



Ml 



M 2 p 



exp 



^ 2 F 2 



A^{s) + -A K (s) 



(9) 



In order to account for the resonance width, here one should not simply replace 
M 2 — s by M 2 — s — iMT(s) in the effective propagator, since this would dou- 
ble count ?sm[Ap(s)] and analyticity would be violated at 0(p 6 ) in the chiral 
expansion. One could avoid the double counting by shifting the imaginary part 
of the loop functions from the exponential to the propagator |31j , but still ana- 
lyticity would be lost. We follow instead a procedure similar to that proposed 
in Ref. [38] for the Kir form factor, in which unitarity and analyticity are pre- 
served. As a starting point we consider a form factor in which the loop functions 
are resummed into the denominator, 



M 2 p 



M 2 p 



M 2 p 



1 + (ArOO + \A K { S j) 



M 2 p 



l+g^ w ^e(A 7T (s) + lA K (s)) -s-iM p T p (s) 



, (io) 



where we have defined the imaginary part of the denominator as —iM p T p (s). 
The energy dependent width is thus given by 



MpS ^ 
:-sm 



96 7T 2 F 2 



Ak(s) + -A K (s) 



(11) 



and from Eq. ^ one has 



sM p 
96 7T F 2 



(s-Ami) aUs) + -9(s~Am 2 K ) <( S ) 



(12) 



which is in agreement with the result obtained from a chiral theory with reso- 
nances [IS] if one assumes the relation Gy = Fy/2. 

The form factor in Eq. (|10p has the correct low-energy behaviour at 0(p 4 ) 
[3T] and leading 0(p 6 ) contributions in xPT [49], and vanishes at short distances 
as expected from the asymptotic behaviour ruled by QCD. As stated, the loop 
functions Ap(s) contain the logarithmic corrections induced by final state in- 
teractions. Now we take into account the fact that the two pion vector form 
factor is an analytic function in the complex plane, except for a cut along the 
positive real axis starting at the threshold for two pion production, s t hr = 4m 2 , 
where its imaginary part develops a discontinuity. From unitarity it can be 
shown [3TJ that the form factor satisfies an n-subtracted dispersion relation 
that involves the scattering phase in the elastic region, for which experimental 
data are available. In the case of n subtractions at s = 0, the dispersion relation 
admits the well-known Omnes solution 



{CXJ 
£l [ ds ' ^ifl 



ie) 



(13) 
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where 

n-l 



logP n ( S ) = ^ a k - (14) 



fc=0 



is the corresponding subtraction polynomial, and Sl(s) is the I = 1, J = 1 two- 
pseudoscalar scattering phase shift. The subtraction constants a k are given 

by 

<*k = -TT log 



ds fc 



(15) 



At least one subtraction is required in Eq. (I13[) to achieve convergence. Here 
the first subtraction constant has been fixed from the normalization Fy(0) = 1, 
which holds with good approximation in view of the conservation of the vector 
current in the isospin symmetry limit. On the other hand, in order to determine 
the scattering phase shift £[(s) [to be used as input in Eq. (fT3")) ] we follow the 
approach in Ref. [55]) taking 

tanffOO = TZJ } , (16) 



where Fy (s) is given by Eq. (ITU1) . In this way, the form factor in Eq. (fTU)) 
trivially satisfies the Omnes relation (|13D for n = 1 and -Fy(O) = 1. This 
form factor should be adequate to reproduce the experimental observations in 
the low energy limit, since by construction it matches xPT results. However, 
beyond this limit one would not expect a sufficiently accurate description of 
the data. Fortunately, the analyticity properties of Fy(s) allow to increase 
accuracy by considering more subtractions in Eq. flT3"l) : for higher n, the weight 
of the dispersive integral at large energies gets reduced, and the corresponding 
information is translated to the subtraction constants [32], which can be taken 
as unknown parameters. In addition, some approach has to be used to deal with 
the phase shift beyond the inelastic two-kaon threshold, where the contribution 
of the dispersive integral is in general still relevant and Eq. (fT3"]) is no longer valid 
(in fact, this happens already at the four-pion threshold, but higher multiplicity 
intermediate states are expected to be phase space suppressed). The goal is 
to obtain a form factor that leads to a satisfactory description of the available 
data, considering just a few subtractions and a phenomenologically adequate 
elastic phase shift. 

On the basis of the previous discussion, we have carried out fits of Fy(s) to 
Belle data from r decays. We find that a good description of the data can be 
obtained with n — 3 subtractions, i.e. taking 



F v {s) = exp 



oil s + — s 



a, „ 2 , t r dsl 6{(s') 



2 nj stht (s>) 3 (s> 



(17) 



In this form factor we have four parameters, namely the subtraction constants 
and the parameters M p and F v that determine the phase shift S\ according 
to Eqs. (fTD|) and (fT5]l . In addition, in order to deal with the phase shift in the 
large energy region, we have distinguished two intermediate (squared) energies 
si and S2- The former is defined as the limit up to which we consider Eq. (|16[) 
to be a reliable description of the phase shift. As stated, we expect this value to 
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be of the order of the inelastic two-kaon threshold, s± ~ 4M|-, or alternatively 
we can consider the limit s± ~ (M p + T p ) 2 , where the effect of the p resonance 
should dominate. In any case one expects si to be about 1 GeV 2 . As we show 
below, this will be supported a posteriori by the good agreement between our 
predictions and the experimental data for 5\ quoted some time ago in Refs. [50j 
Beyond Si, the dispersive integral in Eq. ([TTf should be affected not 
only by inelastic contributions but also by the presence of excited resonance 
states. The other point, S2, indicates the energy at which we assume that the 
phase shift saturates its asymptotic value 5l(s —} 00) = ir, corresponding to the 
existence of a single narrow resonance (5TJ [521 53j3- Here we take S2 — M 2 . In 
the intermediate region, s\ < s < S2, we assume for simplicity a linear behaviour 
of 5\ with s. In order to take into account the uncertainties arising from these 
assumptions, when performing our fits we have considered possible variations of 
the values of s\ and S2, and of the upper integration limit Soo, which is usually 
taken to be in the range [2.25, 00] GeV 2 [32, 38, 54j. The corresponding effects 
on our results have been taken as part of the systematic error of our theoretical 
approach. 

Another aspect to be taken into account is the effect of isospin-violating 
corrections to the pion vector form factor. In general one has to distinguish 
between the neutral and charged pion vector form factors, the latter being de- 
fined by Eq. ([1]) . The corrections can be expanded in powers of the quark mass 
difference and the electromagnetic coupling, in addition to the chiral counting. 
At the leading order, the spectral function for the decay r~ — > ir~ir°v T can be 
written as [55] 

dT(r- -> QV) = ^„, v «(,_JL\ a (,_*L 



d S ~ 384*3 WlKa, ^ M2 j y ^ 

A 3/2 (l,^,^f ) \f + (s)\ 2 G EM (s), (18) 

where Sew = 1 + ( Q: / 7r ) 1°e{Mz/M%) includes the dominant short-distance 
electroweak corrections, and the factor Gem(s) arises from the contribution 
of electromagnetic loops. In the isospin limit one has Sew = Gem = 1, 
A 1 / 2 (l, m^/s, m 2 /s) = cr w (s), and the form factor in Eq. (IT5)l reduces to the 
pion vector form factor, f+(s) — Fy{s). The dominant isospin breaking effect 
in Fy(s) is that arising from phase space, i.e., from considering different masses 
for the charged and neutral pions and kaons in the loop functions. Thus one 
has to replace the functions A 7r (s) and Ak(s) by ^-.^(s) and A K - K o(s), re- 
spectively. Explicit expressions for these functions are given in Appendix A. In 
this way, following the same steps that lead to Eq. (fT7|) we can obtain a dis- 
persion relation for the charged pion vector form factor, Fy + (s). In addition, 
in /+(s) one has to take into account a local electromagnetic correction f^cli' 
which contributes as an additional term in the decay amplitude [55] . One has 
then 

f + (s) = F v + (s) + . (19) 



lr The asymptotic limit of the phase shift obtained from Eq. II16H slightly deviates from this 
value, owing to the linear growth of T p with s. 
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This local electromagnetic correction is given by |55j 

where the scale dependence in the last term cancels that in the logarithms. 
At the scale /i = M p , X(/x) is estimated to be between —2.5 and 4.5 |55j . 
Finally, the loop correction Gem{s) has been computed (including resonance 
contributions) in Refs. [S3] and [57]. We notice that this correction has not been 
taken into account in the extraction of the form factor carried out by the Belle 
Collaboration in Refs. [21 158], 

In order to incorporate the effect of isospin corrections and evaluate its 
significance, we have carried out our fits for the tau decay data considering 
three different situations: 

• I) The limit of exact isospin symmetry, in which /+(s) = Fy(s), where 
the form factor is given by Eq. (|I7I) . 

• II) The inclusion of isospin-breaking corrections at the level of kinematics, 
i.e. considering different masses for the charged and neutral particles in 
the loop functions and the kinematical factors in Eq. (fT5|) . 

• III) The inclusion of all lowest order isospin-breaking corrections, as in 
Eqs. CUB 

As a general result, it is found that we obtain a good fit to Belle data [2] for 
s < 1.5 GeV 2 . Our fits have been carried out with the MINUIT package taking 
the first 30 points (s max = 1.525 GeV 2 , with a bin width of 0.05 GeV 2 ). The 
results are shown in Table [T] First of all, it is worth to notice that — even with 
just a few input parameters — the theoretical curve is able to fit the very precise 
set of experimental data with a x 2 /do/ value close to unity. Hence, it is seen 
that within this energy range the data can be described without the inclusion of 
higher resonant states in the theoretical scheme. On the other hand, it is found 
that the effect of isospin-breaking corrections on the fitted parameters is below 
the 2% level. 





Fit value (I) 


Fit value (II) 


Fit value (III) 


M p [GeV] 
[GeV] 
ai [GeV~ 2 ] 
a 2 [GeV" 4 ] 


0.8431(5)(17) 
0.0901(2)(5) 
1.87(1)(3) 
4.29(1)(7) 


0.8280(4)(14) 
0.0902(2)(4) 
1.84(1)(3) 
4.34(1)(7) 


0.8276(4)(21) 
0.0906(2)(4) 
1.81(1)(2) 
4.40(1)(6) 


x 2 /dof 


1.37 


1.39 


1.56 


T P (M 2 ) [GeV] 


0.207(1) (3) 


0.194(1)(3) 


0.192(1)(4) 



Table 1: Results of our fits. The first and second numbers in brackets correspond to 
the statistic and theoretical systematic errors, respectively. T p (Mp) is obtained using 
the fitted values of M p and F n and is given only for reference. 

Regarding the errors in the fitted parameters, we have quoted separately 
those arising from the fit and the systematic errors coming from the theoretical 

2 For the factor Gem( s ) we have considered the analysis in Ref. I57| . 
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approach. The latter are basically due to the uncertainties in the energy range 
to be fitted, the number of subtractions considered, and the values of s\, S2 
and Soo in the dispersive integral. In order to have an estimation of the effect of 
these uncertainties we have considered the fits for s max in the range [1.325, 1.525] 
GeV 2 , 2 to 4 subtractions, and si, S2, Soo in the ranges [0.95, l.ljj, [M 2 ,oo] and 
[2.25, oo] GeV 2 , respectively. The corresponding results have been quoted in the 
second brackets in Table [TJ while the numbers in the first brackets stand for the 
statistical errors arising from the fit. It is found that M p and F v appear to be 
anticorrelated, and the same happens with the parameters ot\ and ol-2.. 

From the table it is seen that the central values of F n obtained from the fit 
are about two percent below the value of 92.2 MeV quoted by the PDG [59] . The 
difference can be attributed to further theoretical uncertainties, mainly arising 
from the effect of higher order terms in the large Nq expansion. This includes 
corrections to the relations Gy — Fy/2 and GyFy = F%, which have been 
used for the matching between the form factors obtained within the low energy 
xPT theory and the chiral theory with resonances. In fact, as already pointed 
out in Refs. [3J3 [35], the energy-dependent width given by the imaginary part 
of the loop function with F n = 92.2 MeV is not adequately normalized so as 
to reproduce both the experimental data on tttt and Kir tau decay channels. 
Concerning the properties of the p(770) resonance, in order to obtain the corre- 
sponding physical mass and width one should compute the position of the pole 
of the pion vector form factor in the complex s plane, say s po i e . One has 

= K° lC ~ ■ ( 21 ) 

Unfortunately, s po i e cannot be obtained directly from the expression for the 
pion vector form factor in Eq. (1171) . since in general the complex variable s in 
the dispersion relation is not in the same Riemann sheet in which the pole is 
located. We can only obtain an approximate result for the value of s po io by 
looking for the complex pole of the form factor Fy i '- \s) in Eq. (fTTTjl . where the 
values of the parameters M p and F„ are those arising from the fits. Taking 
into account the prescriptions in Refs. [501 EB EH] to deal with the cuts in the 
complex functions entering Eq. (|10[) . we obtain 

M poic = ( 751 3 ± Q g) Me y ^ ppoio = ( 155 6 ± Q g) MeV ( Fit j) . 
M poic = (744.1 ± o.7) MeV , rP° lc = (149.9 ± 0.7) MeV (Fit III) . (22) 

We recall that the form factor Fy^°\s) is used to obtain the input phase 5\(s) 
entering the dispersive integral in Eq. (|17[) . for s < Am 2 K . Our use of the Omnes 
relation is limited by our lack of knowledge of the high energy region, hence 
there can be some difference between the values of M^ olc and rP olc in Eqs. (|2"2")l 
and those corresponding to the actual poles of the form factor. Although the 
shift is expected to be relatively small (the subtractions reduce the weight of 
the high energy part of the dispersive integral), it is difficult to specify the 
corresponding systematic errors, and the values in Eqs. (|22[) have to be taken 
with some care. As a check of consistency, we can compare the phase of the form 
factor obtained from our fit with present experimental data on 8l(s) from Ochs 

3 Best fits are obtained in all cases for s\ ~ 0.98 GeV 2 , in agreement with theoretical 
expectations. 
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et al. [50] , Estabrooks and Martin [53] and Protopopescu et al. [M] . The results 
are shown in Fig.[TJ We quote the data from threshold up to si ~ 1 GeV 2 , i.e. the 
region in which the phase of Fy (s) has been used as input for the dispersive 
integral. In general it is seen that the agreement is very good. It is remarkable, 
however, that our predictions are somewhat below the data in the region of 
very low energies. On the other hand, as shown in the hgure, our results in that 
region are in very good agreement with those recently obtained in Ref. |65) using 
once-subtracted Roy-like equations |66j . Within errors, we also find agreement 
with the values obtained supplementing Roy equations with chiral symmetry 
constraints |67l 168) . This discrepancy between theory and experiment in the 
very low energy region would deserve further tests, taking into account that the 
experimental information corresponds to rather old measurements. 




Vs (GeV) 

Figure 1: Two-pion phase shift Si as function of the mr invariant mass squared. Our 
theoretical expression (red curve) is shown to be in good agreement with experimental 
data (from Ochs et al. [50], Estabrooks and Martin |63] in the s and t channels and 
Protopopescu et al. [64]) up to the opening of the two-kaon threshold, Si ~ 1 GeV 2 . 
At very low energies, where no data are available, our prediction agrees with the results 
of Gara'a-Martfn et al. |65| . 



For comparison, we have also analyzed the results for the pole mass and 
width of the p meson corresponding to the parametrization proposed time ago by 
Gounaris and Sakurai (GS) [BJ, which has been used in the fits carried out by the 
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Belle Collaboration 750. The results obtained by Belle using the normalization 
F£(0) = 1 yield the parameter values M p = (774.9 ± 0.6) MeV, T p = (148.6 ± 
1.8) MeV, which correspond to 

M poie = (761.1 ±0.6) MeV , T p p olc = (142.7 ± 1.6) MeV . (23) 

As an educated guess, we may consider the difference between the results in 
Eqs. (|22p and (j23|) as a conservative estimate of the systematic error in our 
determination of the pole parameters. 

On the other hand, notice that the values for the p mass and width quoted 
by the PDG are 775.49 ± 0.34 MeV and 149.1 ± 0.8 MeV, respectively [59] , 
In fact, these values correspond to the parameters appearing in phenomeno- 
logical amplitudes where the resonances are introduced through BW functions, 
hence they are strongly model dependent Alternatively one can take as the 
relevant resonance properties the above mentioned pole mass and width, which 
should be less dependent on the parametrization. For comparison, in Table [2] 
we show other determinations of the pole mass and width quoted in the liter- 
ature. Though within our approximations we obtain a p pole mass somewhat 
lower than these previous results (presumably owing to the above mentioned 
theoretical uncertainty), a global agreement is found within a two percent level. 
In turn, the values of MP ole and rP olc that we have obtained from the GS fit, 
where the pole parameters can be calculated directly from the corresponding 
phenomenological form factor, are in good agreement with the average results 
from Table H 



Reference 


M polo 


ppolc 


Sanz-Cillero et al. [34] 


764.1 ± 2. 7±f;g 


148.2 ± 1.9±^ 


Ananthanarayan et al. [69] 


762.5 ±2 


142 ±7 


Feuillat et al. [70] 


758.3 ±5.4 


145.1 ±6.3 


Pelaez [71] 


754 ± 18 


148 ± 20 


Zhou et al. [72] 


763.0 ±0.2 


139.0 ±0.5 


Masjuan et al. [73] 


763.7 ± 1.2 


144 ±3 


Results from our fit I 


751.3 ±0.9 


155.6 ±0.8 


Results from our fit III 


744.1 ±0.7 


149.9 ±0.7 


Results from GS parametrization 


761.1 ±0.6 


142.7 ± 1.6 



Table 2: Comparison of different results for the pole mass and width of the p(770) 
meson (values are in MeV). 



Another result that can be derived from our analysis is the value of the point 
satisfying 5\(s^/ 2 ) = 7r /2- This point is often used for the definition of 
the so-called peak mass or visible mass (since it can be directly extracted from 
experimental data). According to our fits, this value is rather stable with respect 
to changes in the above mentioned systematic sources of theoretical errors. We 
find 

v /s^ = (775.0 ±0.2) MeV . (24) 

4 Details of this parametrization are given in the next section, see Eq. 1130 II and below. 
5 This has been pointed out in Refs. 38 54 for the case of the if* (892) resonance, analyzed 
in the context of t~ — > (Ktt)~i> t decays. 
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This is in agreement with the result of Ref. [69 , ^fs^Tpi = (774 ± 3) MeV, 
within a Roy equation analysis of the different partial waves in itit scattering. 
The peak mass has to coincide with the pole mass [62] when the corresponding 
resonance is sufficiently narrow, it lies far from other resonances, and there is a 
negligible background (non-resonant) contribution. 

Finally, we point out that our results for the subtraction constants a,\ and 
Q.i are compatible with previous determinations, although the errors are smaller 
owing to the high quality of present Belle data. We come back to this issue in 
Section |U 



3 Intermediate energy description of Fy(s) 

The approach presented in Sect. [5] has been used to obtain a theoretical de- 
scription of the r — > -khVt spectral function for mr invariant masses up to 
s max ~ 1.5 GeV 2 . Above these energies this description is not adequate, in 
fact, the experimental data are compatible with the presence of excited reso- 
nances. In order to complement the dispersive representation for the pion vector 
form factor proposed in the previous Section, we propose for the energy region 
•Smax ^$ s < M 2 an effective form factor that includes two excited states, namely 
the p'(1465) and //'(1700) resonances. We stress that the dispersive represen- 
tation is insensitive to the modelling of this "intermediate" energy region. The 
dynamics related with the excited resonance multiplets is poorly known, there- 
fore there are no grounds to go beyond a phenomenological parametrization that 
includes several unknown parameters. We extend the form factor in Eq. (|10[) 
as [H] 



M 2 + (a'e 1 *' +a"e i *")s 



Ml 



1 + 96^pT (Ar(s) + \A K (*)) 



M} [l + sCsM*)] 



M 2 p „ [l + sCpt-Ms)} 



(25) 



where the constants C p < and C p " are given by 



C R 



irM R al(M 



(26) 



The resonance masses M p / jP " and on-shell widths T p ^ p n are free parameters of 
this effective form factor. By construction, the off-shell widths of the excited 
resonances behave in a similar way as the p width |43j , considering the two-pion 
states as the dominant absorptive parts of the corresponding self-energies: 



r*(*) 



4(f) 



4rnf 



(27) 



In addition, the form factor includes the coefficients a 1 and a" , which measure 
the relative weight between the contributions of different resonances, and the 
phases <f>' and </>", which account for the corresponding interference. 

Now the unknown parameters can be fitted to Belle data on the r — » mrv T 
spectral function. The quality of the matching between the phenomenological 
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form factor in Eq. (f^5|) and the dispersive representation in Eq. (|10[) can serve as 
a test of the consistency of our approach. The results of our fit for the resonance 
parameters can be translated to the corresponding pole values, leading to 



M poio = ^ 44 ± Q 06 ^ MeV ^ ppoie = (-q 32 ± 0.O6) MeV 
M p ; lc = (1.72 ± 0.03) MeV , T^ 10 = (0.18 ± 0.03) MeV , (28) 



in good agreement with the values quoted by the PDG [SS]. For the coefficients 
and phases we obtain 



a' = 0.08±°;°? 0' = 0.14±°;£° 

a" = 0.03 ± 0.01 <p" = 3.14t°;^° 



(29) 



with a x 2 /dof — 0.26, starting from s = s max . The matching to the form 
factor in Eq. dTUJ) happens at s ~ 1.35 GcV 2 for s max <G [1.3,1.55] GeV 2 . Our 
final curve for the pion vector form factor covering the full range of values from 
threshold to M 2 is shown in Fig. [2] (solid line). The quality of the fits is reflected 
in the good agreement between our results and the experimental data obtained 
by Belle, in particular in the low energy region, where the latter are very precise. 
In addition, it can be seen that the matching at s max ~ 1.35 GeV 2 is smooth, 
which supports the consistency of the phenomenological description proposed 
for the intermediate energy region. In order to appreciate the agreement with 
data with more detail, two close-ups of Fig. [2j corresponding to the low energy 
and the peak regions, are shown in Fig. [3] 

It is worth to point out that the phenomenological form factor in Eq. (|2"5)) 
is qualitatively similar to the GS parametrization [6 mentioned in the previ- 
ous section. Indeed, the GS form factor is built as a sum of Breit-Wigner-like 
functions that keep a nontrivial real contribution in the corresponding denomi- 
nators: 

Fy iGS) (s) = [BW° S { S ) +PBWf, s { S ) +1 BWf,?{ S )} , (30) 

where _ 

M 2 R (1 + d R T R (s)/^) 

BW R (s) = —2 ' — , (31) 

(M R - s) + f R (s) - iM R T r (s) 

and the coefficients /3 and 7 are complex numbers. The energy-dependent widths 
Fij(s) are given, as in our approach, by Eq. (|27[) . while the expression for the 
(real) functions f R (s) can be found in Ref. |B]. The constants d R are chosen 
so that BW R s (0) = 1. As stated, this phenomenological parametrization has 
been used in the fits carried out by the Belle Collaboration [7S], allowing a 
quite successful description of the data throughout the full spectrum. It is 
represented by the dashed curve in Fig. [2] (in the close-ups in Fig. [3] our curve 
and the GS curve overlap, and little differences can only be appreciated in the 
peak region, where our curve shows a slightly better agreement with the data). 
For comparison we also include in Figs. [5] and [3] the result obtained in Refs. [3T] 
and [35] . The latter corresponds to a dispersive representation of the form factor 
in the isospin limit, without the inclusion of excited resonant states (we have 
refitted the parameters according to present Belle data). 
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12 3 



s (GeV 2 ) 

Figure 2: Pion vector form factor Fy(s) compared to Belle data 2. (black dots). 
Solid and dashed lines correspond to our description and the GS parametrization, 
respectively. The dashed-dotted curve stands for the result from Ref. [3T] (for M p = 
775 MeV), while the dotted line corresponds to the dispersive representation in Ref. [32] 
(for qi = 1.83 GeV -2 , a 2 = 4.32 GeV -4 and M p = 774.2 MeV). 




I , I , I , I I I I i I i L 

0.1 0,2 03 0,4 0,4 0,6 0.S 

s (GeV 2 ) s (GeV") 



Figure 3: Two close-ups of Fig. [2] are displayed, corresponding to the low-energy 
region (left panel) and the peak region (right panel). 

4 Low-energy observables 

On the basis of the theoretical approach presented in Sect. [5] we can obtain the 
values of chiral low energy observables. If the expansion of the pion vector form 
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factor in powers of s is parametrized as 

F v (s) = 1 + ~ (r% s + c v s 2 + d v s 3 + ... , (32) 

from Eq. (|17p one has 

(r 2 Y v = Qa 1 , c v = ^(a 2 +a 2 ) . (33) 

Taking into account the results of our fit (case III, i.e. including isospin-breaking 
corrections), we obtain 

(r 2 ) y = 10.86 ±0.14 GeV" 2 , c v = 3.84 ± 0.03 GeV~ 4 . (34) 

These values are indeed in good agreement with almost all previous determina- 
tions made by several authors within various chiral models, see Refs. [2"TI 12"51 1321 
HQ E3 HQ [72 IBS HD IE2]- In order to go beyond the s 2 term in the expansion, 
one can make use of the general relation 

ak = nL , ' (35) 

" J Am \ S 

which allows to determine the subsequent subtraction constants in the Omnes 
expression ()13[) . In this way we obtain a 3 = 29.2 ± 0.2 GeV~ 6 , which leads to 

d v = - (a 3 + 3«ia 2 + al) = 9.84 ± 0.05 GeV" 6 . (36) 

Previous evaluations for this observable have been carried out in Refs. [73] and 
[ST] , leading to 9.70 ± 0.40 GeV~ 6 and 10.18 ± 0.27 GeV" 6 , respectively. In 
order to check the consistency of our procedure we have also calculated the 
constant a>2 from the general relation (1331) . obtaining a 2 = (3.7 ± 0.2) GeV -4 , 
in reasonable agreement with the result of the fit. Notice that, even if Eq. ((35]) 
is exact, one can expect some deviation from the fitted value of a 2 owing to the 
ad-hoc treatment of the phase shift above the inelastic threshold in our analysis 
of the form factor. In the case of a± we cannot trust the result from Eq. ([33)1 
since the slow convergence of the integral provides a large weight to this high 
energy contribution. 



Reference 


r& (M p ) x 10 3 


r&(M p ) x 10 4 


VMD [27] 
0{ P 6 ) xPT [27] 
Pich and Portoles [35] 
Our result 


-0.25 
-0.68(26) 
-0.79(19) 
-0.91(16) 


2.6 

1.50(44) 
1.46(3) 
1.49(1) 



Table 3: Counterterm combinations extracted from (r 2 )y and c v in 0(p 6 ) x?T. 

In addition, the observables (r 2 ) v and c v can be related to two counterterm 
combinations in the 0(p 6 ) chiral Lagrangian, namely r v {M p ) and r V2 (M p ) [2T] . 
which are dominated by the vector resonance contributions (M p ) and ry 2 (M p ). 
Our results for these quantities are quoted in Table |3J showing a good agree- 
ment with the values previously obtained in Ref. [35] and in the 0(p 6 ) xPT 
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fit in Ref. [37]. Within VMD these counterterms can be determined by inte- 
grating out vector resonances in the framework of a chiral effective theory |19j . 
Considering just the contribution of the p meson resonance, within the Proca 
formalism one gets [27] 

r^=^V2^f x f V) rl=^g v f v , (37) 

where fv, gv & n d f x are effective couplings in the chiral Lagrangian with res- 
onances. Our results for rK and ry 2 would lead then to the ratio f x /gv = 
— 2.1 ± 0.5, far from the phenomenological value f x /gv — —0.33 [27]. This indi- 
cates that the role of heavier resonances is crucial in order to describe the 0(p 6 ) 
vector driven contributions in xPT, in agreement with Ref. [32j . 



5 Conclusions 

The high quality data on the pion vector form factor obtained at flavour factories 
demands a correspondingly improved analysis from the theoretical side. In 
order to describe these data keeping the connection with the underlying strong 
interaction dynamics, one can take profit of QCD symmetries to reproduce the 
data in the very low energy domain, and make use of general properties of 
quantum field theory to extend the analysis to higher energies. In this spirit, 
we have presented a dispersive representation of the charged pion vector form 
factor that fulfills the constraints imposed by analyticity and unitarity, and 
reduces to the result obtained within x?T at low energies. 

Our construction is based on the dispersion relation between the form factor 
and the 5\(s) phase shift of elastic tttt scattering. The phase shift is obtained 
from the leading contribution arising in the large- Nc expansion including p(770) 
exchange up the onset of inelasticities, with the further assumption of a smooth 
growth up to the asymptotic value. In this way we obtain a theoretical expres- 
sion for the form factor in terms of four parameters, namely M p , F^, and two 
subtraction constants a\ and ai- The values of these parameters have been 
determined by performing a fit to the very precise Belle data on the r — > irirv T 
spectral function up to a squared 7T7t invariant mass s max — 1.5 GeV 2 , leading 
to the results quoted in Table [1] It is seen that the effect of isospin corrections 
on the parameters lies below the two percent level. From these results we have 
determined the pole values of the p mass and width and the so-called visible or 
peak p mass. We have also obtained the values of low energy observables and 
compared the results with those arising from chiral effective theories. 

In addition, we have addressed the energy region s > s max ~ 1.5 GeV 2 , in 
which the inclusion of excited states is necessary to get a proper description of 
r — > -ktiVt data. For this region we have proposed a phenomenological expression 
for the form factor that takes into account the presence of the resonances p' and 
p" , assuming that the effective propagators behave in a similar way as that of 
the p meson. This allows a good fit to the data, leading to values for the p' and 
p" masses and widths similar to those quoted in previous works. It is seen that 
the curves for the form factor obtained for both energy regions match smoothly 
at s ^ s max . 

As a conclusion, we have seen that the results at low energies supple- 
mented with the leading contributions in the large- Nc expansion are able to 
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provide the input to a dispersive representation of the pion vector form factor 
which fulfills analyticity and unitarity. On this basis, complemented with a phe- 
nomenological description in the high energy region, we have shown that it is 
possible to reproduce the very precise data on t~ — » 7r~7r°^ T decays throughout 
all the phase space. This can be used as input to the new hadronic currents of 
the TAUOLA Monte Carlo generator. On the other hand, our fits lead to the 
parametrization of the charged pion vector form factor, thus the comparison 
with a precise determination of the neutral form factor would provide robust 
information on the mr contribution to the muon anomalous magnetic moment. 
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Appendix 

The explicit form of the loop functions Apq(s) can be obtained from Ref. |47j . 
One has 

= _ 192 ir* [B Mp q Q0-W)] , (38) 
s 

where Mpq(s) and Lpq(s) can be written in terms of new functions Epq, Apq, 
kpQ, Jpq and J PQ as 

1 - Ap Q .1 i 

MpQ ^ = 12~^ S ~ 2EpQ ) Jp ^ s ) + Y^ JpQ ^'~ Q k PQ + 
A 2 

Lpq{s) = -^Jp Q ( s ). (39) 

The new functions Epq and Apq are defined by £pg = mp + nig, Apq = 
mp — rriQ, while fcpQ includes the renormalization scale /x: 

k PQ = ^MMP - MQ) , (40) 

ApQ 



rri r, , / m 



where 

^ = £hi>'°z{w) (41) 

(we have taken \i = M p , as in the isospin symmetric case). Finally, the functions 
Jpq and Jpq are given by 

Jpq(s) = Jpq(s)- sJ pq (0) 
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Jpq{s) 



32 tt 2 



^PQ ^PQ \ , ( 



log 



V ( (S+V) 2 -&PQ 



(42) 



where ^ = A 1 ^ 2 (s, mf,, mg). We note finally that 



» J ^°)^(Ag + 2 ^f l0g sfj ■ <43) 
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